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ABSTRACT 

By making use of Herschel-PEP observations of the COSMOS and Extended Groth 
(Z~) Strip fields, we have estimated the dependence of the clustering properties of FIR- 

^__l selected sources on their lOO^fm fluxes. Our analysis shows a tendency for the clustering 

y—i strength to decrease with limiting fluxes. By assuming a power-law slope with 7 = 1.8 

£ — . for the two-point correlation function £(r) = (r/ro)~ 7 , we find: ro(iSioo^m z? 8 mJy) = 

r—4- 4.3^q7 Mpc and ^(.Sioo^m ^ 5 mJy) = 5.8^ 2 'o- These values convert into minimum 

(~^> halo masses M m i n ~ 10 116 M for sources brighter than 8 mJy and M m ; n ~ 10 124 

CO Mq for fainter, (Sioo^m ^ 5 m Jy galaxies. We show such an increase of the clustering 

strength to be due to an intervening population of z ~ 2 sources, which are very 

*>. strongly clustered and whose relative contribution, equal to about 10% of the total 

• ^h counts at Sioo^m ^ 2 mJy, rapidly decreases for brighter flux cuts. By removing 

/\ such a contribution, we find that z < 1 FIR galaxies have approximately the same 

5-j clustering properties, irrespective of their flux level. The above results were then used 

to investigate the intrinsic dependence on cosmic epoch of the clustering strength of 
dusty star- forming galaxies between z ~ and z ~ 2.5. This was done by comparing 
our dataset with IRAS in the local universe and with sources selected at 160/im in 
the GOODS-South at z ~ 2. In order to remove any bias in the selection process, the 
adopted sample only includes galaxies observed at the same rest-frame wavelength, A ~ 
60/im, which have comparable luminosities and therefore star-formation rates (SFR> 
100 M Q /yr). Our analysis shows that the same amount of (intense) star forming 
activity takes place in extremely different environments at the different cosmological 
epochs. For z < 1 the hosts of such star forming systems are small, M m j n ~ 10 11 M Q , 
isolated galaxies. High (z ~ 2) redshift star formation instead seems to uniquely take 
place in extremely massive/cluster-like halos, M m ; n <~ 10 135 M Q , which are associated 
with the highest peaks of the density fluctuation field at those epochs. 

Key words: galaxies: evolution - galaxies: statistics - infrared - cosmology: observa- 
tions - cosmology: theory - large-scale structure of the Universe 
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1 INTRODUCTION 

Investigations of the Large Scale Structure (LSS) traced by 
selected classes of extra-galactic sources, combined with in- 
creasingly more refined models for its description, has en- 
abled in the recent years the derivation of important infor- 
mation on some of the physical properties characterizing the 
objects which produce the clustering signal. One of these key 
quantities is the halo mass and its evolution with look back 
time, which for instance allows one to derive cosmological 
connections between different populations of extra-galactic 
sources observed in different wavelengths at different epochs. 

Within this framework, amongst the more interesting 
classes of sources are those which are undergoing some kind 
of active process such as AGN or star formation. Since 
the first high redshift observations of the two Hubble Deep 
Fields (Williams et al. 1996; Casertano et al. 2000), all clus- 
tering analyses performed on optical/near-IR surveys which 
were at the same time wide enough to ensure statistical sig- 
nificance and deep enough to probe high redshift regimes, 
have highlighted a strong evolution of these two classes 
of sources, whereby activity shifts to smaller and smaller 
haloes/galaxies as one moves forward in time. The critical 
redshift for this transition seems to be around z ~ 2 for 
star-forming galaxies (e.g. Magliocchetti & Maddox 1999; 
Arnouts et al. 2002; Foucaud et al. 2010; Hartley et al. 2010) 
and z > 3 for QSOs (e.g. Porciani, Magliocchetti & Norberg 
2004; Shen et al. 2007; Ross et al. 2009; Shen et al. 2009). 
Perhaps not surprisingly, these two redshift values bracket 
the peak of cosmic activity both in terms of AGN and star 
formation (Madau et al. 1996). 

This pattern, referred to as "downsizing" (Cowie et al. 1996; 
Heavens et al. 2004), is also observed in an increasing num- 
ber of independent works (see e.g. Treu et al. 2005; Cimatti, 
Daddi & Renzini 2006, Saracco et al. 2006; Bundy et al. 
2006; McLure et al. 2006 just to mention a few), all adopt- 
ing techniques which differ from clustering analyses. 

The advent of the Spitzer telescope has marked another 
milestone in our knowledge of the population of dusty (and 
therefore optically dim) active galaxies. In particular, for 
the first time Spitzer has detected a population of Ultra- 
Luminous Infrared Galaxies (ULIRGS, L > 10 12 Lq) in the 
redshift range (1.6 < z < 2.5) with extremely high clustering 
lengths (ro ~ 15 Mpc - Farrah et al. 2006; Magliocchetti et 
al. 2007; 2008; Brodwin et al. 2008; Starikova et al. 2012). 
Taken at face value, these correlation lengths in the local 
universe correspond to those of groups-to-poor clusters of 
galaxies (see e.g. Guzzo et al. 2000; Estrada et al. 2009), 
and might be interpreted with a roughly one-to-one corre- 
spondence between the overwhelming majority of ULIRGS 
at z ~ 2 and z < 1 clusters (e.g. Brodwin et al 2007). 

The Herschel telescope allows the scientific commu- 
nity to further extend its knowledge on the formation and 
evolution of dusty sources up to very high redshifts. First 
attempts to measure the clustering properties of Herschel 
galaxies have been presented by Maddox et al. (2010), 
Cooray et al. (2010) and Amblard et al. (2011) for sources 
detected with the SPIRE (Griffin et al. 2010) instrument 
at the wavelengths A = 250, 350 and 500/im. However, the 
above estimates can at most be considered as tentative since 
in the first place they do not agree with each other, and most 
importantly, because the de-projected quantities which are 



reported in some of these works are solely based on models 
for the redshift distribution of SPIRE sources at the various 
wavelengths. 

A step forward has been made possible by the Herschel- 
PEP survey, which provides deep observations taken with 
the PACS instrument (Poglitsch et al. 2010) at the wave- 
lengths A = 70, 100 and 160/im on a number of well stud- 
ied fields such as COSMOS, GOODS-North and South, the 
Lockman region, the Extended Chandra Deep Field and 
the Extended Groth Strip. While probing less sources than 
those collected in the mainly SPIRE-based programs, there 
are undoubt advantages connected with the exploitation of 
PEP data: on the one hand, the shorter wavelength range 
implies a smaller angular resolution and therefore an eas- 
ier search for counterparts of PEP-detected sources. On the 
other hand, observations of well known and studied regions 
such as the two GOODS fields or COSMOS, not only allow 
identification and characterization of FIR galaxies at the 
different wavelengths, but also provide spectroscopic and/or 
phometric redshift determinations for the overwhelming ma- 
jority of such sources. 

By making use of PEP observations of the GOODS- 
South field, Magliocchetti et al. (2011) managed for the first 
time to provide direct (i.e. purely based on data) measure- 
ments of the three-dimensional clustering properties of far- 
infrared sources up to z ~ 3. These authors report a cluster- 
ing length ro — 6.5 Mpc for a sample of galaxies with lOO^im 
fluxes brighter than 2 mjy. Possibly more importantly, they 
also find a tremendous increase of the clustering strength for 
galaxies residing at z ~ 2. Indeed, for this population the 
reported values for the correlation length are ro = 17 — 19 
Mpc, which set these galaxies within very massive, cluster- 
scale dark matter halos. 

The results presented in Magliocchetti et al. (2011) are 
in agreement with those found in the literature for actively 
star-forming galaxies selected at redshifts z ~ 2 to z ~ 3 (cfr. 
beginning of this Introduction). However, although some 
steps forward have been taken in order to minimize its ef- 
fect (see e.g. Magliocchetti et al. 2008; Focaud et al. 2010; 
Lin et al. 2012), all the conclusions on a possible intrinsic 
evolution of the clustering strength with look-back time are 
hampered by the well known 'Malmquist bias' issue, since 
in flux limited surveys the more distant sources are also the 
more luminous ones which, at least in the visible part of the 
spectrum, are generally found to be intrinsically more clus- 
tered than fainter galaxies. Furthermore, a survey which is 
carried out at a fixed wavelength/waveband will probe dif- 
ferent portions of a source's electromagnetic spectrum at 
the different redshifts. As a result of this second effect, a 
heterogeneous mixture of galaxy types might end up in a 
survey sample as one moves further back in time, prevent- 
ing to draw any conclusion on the evolution of a specific 
population with redshift. 

The present paper aims at overcoming these two poten- 
tial selection biases so to provide some information on the 
intrinsic evolution of the clustering properties of FIR galax- 
ies. This will be done by estimating the two-point correlation 
function w(8) of PEP-selected galaxies with different fluxes. 
In order to do so, our analysis will include some of the deep 
fields probed by the PEP survey such as COSMOS, the Ex- 
tended Groth Strip and GOODS-South. The work is orga- 
nized as follows: Section 2 introduces the PEP catalogues. 
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Figure 1. Normalized distribution of PEP sources as a function 
of 100//m flux in the Extended Groth Strip (EGS), COSMOS and 
GOODS-South fields. 



Section 3 presents the results for the two-point correlation 
function of PEP sources as a function of their fluxes both in 
projected (§3.1) and in three-dimensional (§3.2) quantities, 
while Section 4 provides some information on the physical 
properties of FIR galaxies such as the halo masses of their 
dark matter hosts. Section 5 uses the results found in §3 
and §4 to investigate the clustering properties of relatively 
local, z < 1, PEP-selected galaxies as a function of their 
flux. This sample will then be compared with others taken 
from the literature for sources selected at the same rest- 
frame wavelength as the ones considered in this work and 
with comparable bolometric luminosities, so to investigate 
the evolution of the properties of dusty star-forming galaxies 
in the redshift range z ~ [0—2.5] in the absence of systematic 
observational biases. Section 6 summarizes our conclusions. 
Throughout this paper we assume a ACDM cosmology with 
H = TOKms^Mpc" 1 , D. = 0.27, Q A = 0.73 and a 8 = 0.8. 



2 THE DATA 

The COSMOS and Extended Groth Srip (EGS) regions have 
been observed by the PACS (Poglitsch et al. 2010) instru- 
ment onboard the Herschel Space Observatory (Pilbratt et 
al. 2010) as a part of the PACS Evolutionary Probe (PEP, 
D. Lutz et al. 2011) Survey, aimed at studying the proper- 
ties and cosmological evolution of the infrared population 
up to redshifts z ~ 3 — 4. We refer to the Lutz et al. (2011) 
and Berta et al. (2010) papers for further information on the 
survey and fields. 

The total number of sources detected at lOO/xni in the 
blind (i.e. without the use of priors, in order to minimize the 
possible bias effects introduced by a double selection) cata- 
logues at the > 3<r confidence level is 5360 in the COSMOS 
region and 1012 in the EGS. The corresponding 100/im flux 
limits are ~ 3 mJy for COSMOS and ~ 2 mJy for EGS. We 



note that in this work we will only make use of observations 
taken at 100/im since, as it was shown by Magliocchetti et al. 
(2011), clustering results obtained at 160/^m were virtually 
identical to those derived at 100/im. 

The distribution of sources as a function of flux in 
these two fields is shown in Figure [T] and it is compared 
with that of 100/^m-selected galaxies observed by PEP in 
the GOODS-South (GOODS-S), to date the region covered 
by the deepest Herschel data. As it is possible to appreci- 
ate from the Figure, the EGS field is > 80% complete for 
fluxes Sioo/jm =? 5 mJy, while COSMOS only reaches the 
above completeness level for Sioo^m ^ 8 mJy. These figures 
are also in agreement with those estimated from simulations 
(cfr. Lutz et al. 2011 and Berta et al. 2010). 

Since clustering studies need high completeness levels, 
in the following analysis we will only concentrate on those 
sources belonging to the blind PEP catalogues which have 
WOfim PACS fluxes above the ~ 80% completeness level. 
Applying the above flux cuts, the number of sources are 
respectively 3913 in COSMOS (Sioo M m > 8 mJy) and 751 in 
EGS (Sioo M m #= 5 mJy). Their projected distributions onto 
the sky are shown in the two panels of Figure [2] 



3 CLUSTERING ANALYSIS 

3.1 The Angular Correlation Function 

The angular two-point correlation function w(9) gives the 
excess probability, with respect to a random Poisson distri- 
bution, of finding two sources in the solid angles SQi SH2 
separated by an angle 9. In practice, w(9) is obtained by 
comparing the observed source distribution with a catalogue 
of randomly distributed objects subject to the same mask 
constraints as the real data. 
We chose to use the estimator (Hamilton 1993) 

DD-RR , 



w(0) = 4 x 



{DRY 



(1) 



where DD, RR and DR are the number of data-data, 
random-random and data-random pairs separated by a dis- 
tance 6. 

Given the requirement for ^ 80% completeness in all 
the considered samples (cfr Section 2) and the considerable 
uniformity of data sampling throughout the fields, we have 
estimated w(9) by simply generating random catalogues of 
about 100K sources so to cover the whole surveyed areas mi- 
nus the outermost regions which showed irregular data cov- 
erage. Within the considered areas, the number of 100/im- 
selected sources is 3005 in the COSMOS field and 647 in the 
EGS. 

w(ff) in eq. |l]) for the above objects was then estimated 
on angular scales ranging from 10 -3 degrees up to values 
of 9 ~ 0.1 — 0.5 degrees, depending on the geometry of the 
two fields since the upper limit cannot be larger than about 
half of the maximum scale probed by a survey. We also note 
that we did not correct for incompleteness as a function of 
IR flux when counting pairs, as its effect would be much 
smaller than the large uncertainties associated with the 
estimates of w(9), due to the relatively small number of 
sources belonging to the considered fields. 
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Figure 2. Projected distribution of sources for the fields considered in our analysis. Left-hand side is for COSMOS, while the right-hand 
side represents the EGS. In order to ensure > 80% completeness, the first field only includes sources brighter than Sioo/jm = 8 mjy, 
while the second one is obtained for a flux-cut SioOfim = 5 mjy. 
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Figure 3. Projected two-point correlation function for 100-/^m-selected sources in COSMOS (Sioo/jm jS 8 mjy, left-hand panel) and 
EGS (SioOfim ^ 5 mjy, right-hand panel). The short-long dashed lines show the best power-law fit, obtained for 7 = 1.8 (fixed) and for 
different values of the amplitude A, while the dashed lines are the best fits to models obtained for a varying minimum halo mass (see 
text and Table 1 for details). 



Figure [3] then shows our results for the angular cor- 
relation function in COSMOS (SiooMm =S 8 mjy) on the 
left-hand panel and EGS (Sioo^m ^ 5 mjy) on the right- 
hand panel. Due to the angular resolution of the Herschel 
mirror at the considered wavelength (~ 8") which would 
cause source blending on lower angular scales, the results 
in the smallest 9 bins of the angular correlation function 



merely represent lower limits and therefore have been indi- 
cated with upward arrows. 

The error-bars have been obtained from jack-knife re- 
sampling. In the first case, the COSMOS field was divided 
into 9 quadrants of approximately the same area and the cor- 
relation function w(6) was calculated for nine realizations of 
the data, obtained each time by omitting one quadrant. Er- 
rors were then obtained from the variance in w. The same 
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Figure 4. Normalized redshift distribution of 100/jm-selected 
sources in the GOODS-South field (Sioojim 1? 2 mjy), for the 
EGS field ( Sioo M m > 5 mjy) and COSMOS (Sioo^m > 8 mjy). 



approach was used on the EGS, where the whole field was 
divided into 10 quadrants. 

If we then assume the standard power-law form for 
w(0) = A6 1 " 1 , we can estimate the parameters A and 7 
by using a least-squares fit to the data. Given the large er- 
rors in w, in all cases we chose to fix 7 to the standard value 
7 = 1.8. Although somewhat arbitrary, this figure and its 
assumed lack of dependence on redshift is justified by Large 
Scale Structure observations of large enough samples of high 
redshift sources so to allow for a direct estimate of the slope 
of w(6) at different look back times (e.g. Porciani, Maglioc- 
chetti & Norberg 2004; Le Fevre et al. 2005). The small areas 
of both the COSMOS and EGS fields introduce a negative 
bias through the integral constraint f w eBt d£lid£l2 = 0. We 
allow for this by fitting to A6 1 ' 1 — C, where C = 1.93 in 
the case of COSMOS and C = 2.54 for EGS were found by 
numerical integration as in Roche & Eales (1999). 

The long-short dashed lines in Figure 13] represent the 
best power-law fits to all the various data-sets. The asso- 



ciated best-fit values for the amplitude are: A(S-_ 



100/jm ;» 

8 mJy) = [1.0 ± 0.3] • 10" 3 for the COSMOS sample and 
A(5ioo M m > 5 mJy)=[l.l+J5;J] ■ 10~ 3 for the EGS sample. 
They are summarized in Table 1. The above values for the 
clustering amplitude are comparable with those obtained 
by Magliocchetti et al. (2011) for the sample of Sioo^m ^ 2 
mJy-PEP selected galaxies in the GOODS-S field. We also 
note that the low clustering signal observed by Maddox et al. 
(2010) for galaxies selected at 250pm within the Herschel- 
ATLAS consortium is compatible with our PEP results. 



3.2 Relation to spatial quantities 

The standard way of relating the angular two-point correla- 
tion function w(6) to the spatial two-point correlation func- 
tion £(r, z) is by means of the relativistic Limber equation 



w(6) 



^ ^ F -2 {x)x ^2 (x) ^ z)dxdu 

[f™ F-i(x)x 2 $(x)dx] 2 



(2) 



where x is the comoving coordinate, F(x) gives the correc- 
tion for curvature, and the selection function Q(x) satisfies 
the relation 



N-- 



$(x)F (x)x dx - 



1 



N{z)dz, 



(3) 



in which N is the mean surface density on a surface of solid 
angle Q, s and N(z) is the number of objects in the given 
survey within the shell (z, z + dz). 

Both COSMOS and EGS have reliable and statisti- 
cally complete catalogues of source redshifts (either spec- 
troscopic or photometric; D. Rosario and O.Ilbert private 
communication). The redshift distributions of PEP-selected 
sources brighter respectively than 5 mjy and 8 mjy in 
these two fields are shown in Figure HI They are compared 
with the N(z) derived for GOODS-S sources brighter than 
SiooMm = 2 mjy presented in Magliocchetti et al. (2011). 
Two main features can be observed in this plot: the first 
one is a systematic trend for brighter sources to be found 
at increasingly lower redshifts (< z >~ 0.65 for sources 
brighter than 8 mjy, < z >~ 0.85 for Sioo^m ^ 5 mjy and 
< z >— 0.92 for 5'ioo/jm #= 2 mjy). The second noticeable 
feature is a progressive disappearance of the 1.7 < z < 2.6 
population with brighter limiting fluxes. We will discuss this 
second point and its consequences on the clustering results 
in more detail in Section 5. 

If we then adopt a spatial correlation function of the 
form £(r,z) — {r/ro)" 1 ' 8 and we consider the redshift 
distributions at the various flux limits as presented in 
Figure [4] from the observed angular correlation function 
and for the adopted cosmology we obtain: r Mm (S^8 
mJy)=4.3±^ Mpc and ^ 00 " m (S> 5 mJy)=5.8J:2:o M P C > 
where both quantities are comoving and the first result refers 
to COSMOS sources, while the second one to the EGS sam- 
ple. We remark that the choice for a fixed value of the quan- 
tity 7 = 1.8, although a sensible one for describing the data 
presented in Figure 3, implies that in the process of data- 
fitting the covariance between the amplitude of the correla- 
tion function and its slope has been ignored. This implies 
that the errors on ro presented above might have been un- 
derestimated. 

From the above results it appears a tendency for 
brighter sources to be less clustered than fainter ones. This 
trend is confirmed if one goes to even fainter fluxes and in- 
cludes the results of Magliocchetti et al. (2011) for Sioo^m ^ 
2 sources in the GOODS-S field (ro°° Mm = 6.3±i'.3 Mpc). 

We also note that the results for the brightest (and 
on average closest, cfr Figure BJ infrared galaxies are vir- 
tually indistinguishable from those obtained locally for 
star-forming objects selected via various methods. Indeed, 
the COSMOS correlation length mirrors those derived for 
optically-selected late-type galaxies (e.g. ro = 5.2 ± 0.4 
Mpc; Madgwick et al. 2003 - r = 4.5 ± 0.2; Zehavi et al. 
2011), for UV-selected galaxies both locally (ro = 4.6 ± 0.6 
Mpc at z~ 0.2; Heinis et al. 2007) and at higher redshifts 
(r = 4.9 ± 0.5 Mpc at z ~ 0.9; Heinis et al. 2007) and 
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Figure 5. Comoving correlation length as a function of 100/xni 
flux as seen in various fields observed by the Herschel-PEP sur- 
vey. The point at 8 mJy is obtained for COSMOS sources, while 
that at 5 mJy is derived for objects belonging to the EGS field. 
The result at 2 mJy is taken from Magliocchetti et al. (2011) for 
sources in the GOODS-S field. 



Figure 6. Minumum dark halo mass as a function of 100/^m flux 
as seen in various fields observed by the Herschel-PEP survey. The 
point at 8 mJy is obtained for COSMOS sources, while that at 5 
mJy is derived for objects belonging to the EGS field. The result 
at 2 mJy is taken from Magliocchetti et al. (2011) for sources in 
the GOODS-S field. 



for IRAS-selected sources (ro = 5.4 ± 0.2 Mpc; Saunders, 
Rowan- Robinson & Laurence 1992). 

Our local results are also in good agreement with the 
very recent estimates of van Kampen et al. (2012), obtained 
for a local (z ^ 0.3) sample of 250/im-selected galaxies taken 
from the Herschel- ATLAS survey so to posses reliable spec- 
troscopic redshifts (ro = 5.62 ± 1.14 Mpc), although their 
result is likely boosted by the presence of a structure found 
at z — 0.164. Indeed the value reported by van Kampen et 
al. (2012) for the correlation length in the lowest redshift 
bin reads r = 3.29 ± 1.10 Mpc. 

A visual representation of the dependence of the clus- 
tering length ro on 100/im flux is provided in Figure [5] A 
more detailed discussion of its implications will be presented 
in Section 5. 



4 CONNECTION WITH PHYSICAL 
PROPERTIES: THE HALO BIAS 
APPROACH 

The easiest way to relate the clustering trend of a popula- 
tion of extra-galactic sources with their physical properties 
is by means of the Halo Bias approach (Mo & White 1996; 
Sheth & Tormen 1999) which describes the observed clus- 
tering signal in terms of the product between the two-point 
correlation function of the underlying dark matter distribu- 
tion £dm (r, z) and the square of the so-called bias function 
b(Mmi n ,z), which at a given redshift solely depends on the 
minimum mass M m i n of the haloes in which the detected 
astrophysical objects reside via the relation: 



More sophisticated models have been introduced in the 
recent years in order to provide a more realistic descrip- 
tion of the observed clustering properties of extra-galactic 
sources. One of the most successful ones, the so-called Halo 
Occupation Model (HOM, e.g. Scoccimarro et al. 2001) re- 
lates the clustering properties of a chosen population of 
galaxies to the way such objects populate their dark mat- 
ter haloes, i.e. relaxes the somehow unrealistic assumption, 
implicit in the halo bias model, of having a one-to-one cor- 
respondence between a galaxy and its dark matter halo. 
However desirable, such an analysis is unfortunately not pos- 
sible on our datasets due to a relatively poor signal-to-noise 
ratio. In fact, the HOM is a three parameter model: to- 
gether with the minimum mass of the dark matter halo it 
also depends on the the quantities a and iVo according to 
the expression 



JV(Af) =N (M/M a 



if M > M n 



£(r, z) = ^dm{r, z) ■ b(Mm in , z) 



(4) 



where N(M) is the number of galaxies within a halo of some 
mass M. And a three-parameter fit will only be possible 
when larger fields will be included in the analysis. For the 
time being, we notice that the quantities M m i n and a are 
covariant, so that within the HOM framework, higher val- 
ues for q in general correspond to lower values for M m i n . 
This implies that in the presence of multiple halo occupancy, 
the values for M m j n found within the HOM scenario will be 
(slightly) lower than those obtained via the halo bias model 
adopted in this work. 

Having made this necessary digression, we can then 
go back to the Halo Model framework and compute the 
theoretical angular two-point correlation function w(8) t h 
predicted by this model starting from eq.Q and projecting 
it once again by following eq.pj with the iV(z)'s provided 
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Table 1. Properties of the fields observed by the PEP survey at 100 /im and included in our analysis: a) field name; b) limiting flux 
for ^ 80% completeness; c) total source number above 80% completeness; d) amplitude A of the angular two-point correlation function 
w(8); e) comoving correlation length ro; f) minimum dark halo mass M m j n . 



Field Name 



< ' tin ti 

O 100/Jm 



[mjy] Number sources 



r [Mpc] Logio M min [M ] 



COSMOS 
EGS 



3005 
647 



(i.oig;|) x io- 



(1.1+^) x 10- 



4 3+0-7 
t>.«_ 2 . 



11.6 
12.4 4 



+0.4 
0.7 



11.6i ^ and 



in Figure [4] 

The best values for M m i n at the two limiting fluxes of 5 
mjy and 8 mjy were then found by a x 2_n t to the observed 
clustering signal. This gives: 
Logi [M min /M Q ](S ^ 8 mJy) = 

Logio[M m in/M Q ](S' ^ 5 mJy)=12.4lo;g, where again 
the first result refers to the COSMOS field, while the sec- 
ond one to the EGS sample. These results are summarized 
in Table 1. The corresponding best-fit w(6) t h are shown in 
Figure [3] by the dashed curves. 

The values obtained in our analysis, especially those de- 
rived for the brightest and, as already stressed, on average 
closest sources, are again in excellent agreement with results 
found in the literature for the typical masses of halos locally 
hosting star-forming galaxies. Just to mention a few, Heines 
et al. (2007) quote bias levels which, by using the Sheth & 
Tormen (1999) bias model, convert to minimum halo masses 
of Log lo [M min /M ] ~ 10.7 and Log lo [Af min /M ] ~ 11.5, re- 
spectively for their low-redshift and z ~ 0.9 samples. And 
while the above figures were derived for UV-selected galax- 
ies, very similar results are also obtained by Magliocchetti 
& Porciani (2003) when analysing the clustering properties 
of a sample of star-forming galaxies extracted from the 2dF 
survey (Logio[M m i n /M ] ~ 11) and also by Zehavi et al. 
(2011) for their sample of blue galaxies drawn from the sev- 
enth release of the SDSS (Logio[M min /M ] ~ 11). 

Furthermore, although we have not made use of the 
HOM model to estimate our best-fit parameters, we can 
still rely on its theory to draw some conclusions on the 
sources which produce the clustering signal presented in Fig- 
ure [3] In fact, both the estimated angular correlation func- 
tions, irrespective of the particular field considered in this 
work or of the flux cut, present a flattening (coinciding in 
most cases with negative values of w(0)) on angular scales 
0.004 < 9 < 0.01 degrees. This can be easily explained as a 
further evidence for the well known "avoidance behaviour" 
of local star-forming galaxies which prefer to stand isolated 
and not have companions within the same dark matter halo 
(cfr Magliocchetti & Porciani 2003 or Zehavi et al. 2011 just 
to mention a couple of works). 

As a last point, we show in Figure [6] the dependence 
of the minimum halo mass M m j n on 100/im flux. As in the 
case of the estimated correlation lengths, there is a trend for 
sources brighter than 8 mjy to be hosted in smaller halos. 
This is also further confirmed by the Magliocchetti et al. 
(2011) results which can be used to extend the analysis to 
even fainter, Sioo^m ^ 2, sources. The next Section will 
discuss the implications of both this and other findings of 
our work. 



5 THE REDSHIFT EVOLUTION OF 
CLUSTERING OF FIR GALAXIES 

5.1 Dealing with redshift dependence 

In the previous sections we have found a tendency for the 
clustering signal of PEP sources selected at 100/im to de- 
crease as one moves to higher flux limits. This in turn implies 
that brighter galaxies are hosted by smaller haloes, with val- 
ues ranging from ~ 10 M at the brightest fluxes probed 
by our analysis up to ~ 1O 12 ' 5 M for a ^ 5 mjy flux cut. In 
Section 3.2 we have also learned that, in the case of 100/^in- 
selected galaxies, higher fluxes on average imply more local 
sources (cfr. Figure [4b. Putting together these two pieces of 
information in a straightforward way then would imply a 
clustering strength which increases for increasing redshifts, 
at least up to z ~ 1. But is that so? 

As already extensively discussed, Figure pi clearly shows 
that the contribution of the high, z ~ 2, population to 
the total number of PEP sources, basically negligible for 
Sioojjm ^ 8 m Jyi tends to steeply increase when moving 
to fainter fluxes. Magliocchetti et al. (2011) have shown 
this z ~ 2 population to be a very extreme one, being 
very strongly clustered (ro ~ 19 Mpc) and hosted in very 
massive/cluster-like halos (M > 1O 13,5 M ). So the question 
one might want to ask is: how much of the measured cluster- 
ing signal for fluxes < 8 mjy is solely due to the high-redshift 
population and how much comes from more local sources? 
In other words: can we separate within our data-sets the 
clustering signal stemming from high-z sources from that 
originating from low-redshift ones? 

Indeed, this is possible. In fact, in our case, the observed 
correlation function produced by the two low-z (A) and high- 
z (B) populations can be written as: 



£(r,*) = 






for 
for 



2 < 1.7 
1.7 < \z < 2.6, 



(5) 



so that the observed projected correlation function w{9) will 
be 



){9) = n 2 A WA(0) + n B w B (9) 



(6) 



(see e.g. Magliocchetti et al. 1999), where ua and n B are 
the relative space densities of the two populations at the 
different flux cuts, normalized so that ua + n B = 1. Note 
that, since we want to investigate the clustering behaviour 
of objects located at different redshifts, the cross-correlation 
signal nAn B \wa{9)w b (9)] in equation fej) is zero. 

We take the correlation function £g of lOO^m-selected 
high-z galaxies from Magliocchetti et al. (2011) and estimate 
the values of ua and n B from the redshift distributions pre- 
sented in Section 3.2. Specifically: £ B = (r/ro)~ 18 , with 
ro = l^U; n A (S > 8 mjy) ~ 0.96, n B (S ^ 8 mjy) ~ 0.04; 
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Table 2. Properties of a homogeneous sample of star-forming galaxies observed by the Herschel-PEP and IRAS surveys at A ~ 60^tm 
- rest frame. The first column reports the field name and the considered redshift range. S m i n is the mjy flux cut, ro the correlation 
length and M m i n the logarithm of the minimum halo mass. < z > indicates the mean redshift of each sample, < Arf > the rest-frame 
wavelength, while <L m ; n > and <SFR m ; n > respectively indicate the average minimum bolometric luminosity and the average minimum 
star formation rate. The first five lines are the results of the analysis presented in this work, while the sixth one refers to the sample 
of 160/^m-selected galaxies in the GOODS-S field as obtained by Magliocchetti et al. (2011) and the last one to the results obtained by 
Saunders et al. (1992) for IRAS galaxies. The redshift distribution for IRAS galaxies was taken from Saunders et al. 1990. The superscript 
" 1" illustrates the clustering results obtained by removing the high-z clustering signal by means of equations [5] and [6] while " 2" those 
derived from a direct analysis of the z < 1.7 sample (see text for details). 



Field Name 



r [Mpc] 



M„ 



[M e 



< z > 



< A RF >H <L min >[L ] <SFR min > [M /yr] 



COSMOS[100Atm ] (z < 1.7) 
EGS[100/im ] (z < 1.7) 1 
EGS[100/im ] (z < 1.7) 2 
GOODS-S [100/mi ] (z < 1.7) 1 
GOODS-S [100/im] (z < 1.7) 2 
GOODS-S [160/rni] (1.7 ^ z s£ 2.6) 
IRAS[60Aim] (z ~ 0) 



8 


4 1+0-8 


11 I" 1 " 04 
u - i -0.T 


5 


5 0+ 2 ' 2 
°' U -3.3 


11.9+°-* 


5 


r 1+1.4 


11.8iH 


2 


r l + l.l 
°' i -2.6 


12.0^1 


2 


4 3+1-8 


11 5 +0 ' 9 


5 


1 7 4+ 2 - 8 
1( -*-3.1 


in 7 +0.3 
lci -'-0.4 


600 


5 4+ ' 2 


11 4 +0 - 2 



0.56 ±0.36 


67 ± 14 


0.68 ±0.39 


62 ± 14 


0.68 ±0.39 


62 ± 14 


0.70 ±0.35 


61 ± 13 


0.70 ±0.35 


61 ± 13 


2.1 ±0.3 


52 ±5 


~0.02 


60 



11.6 ±0.6 


72 ± 106 


11.6 ±0.5 


72 ±84 


11.6 ±0.5 


72 ±84 


11.2 ±0.5 


29 ±27 


11.2 ±0.5 


29 ±27 


12.1 ±0.2 


226 ±113 


11.0 ±0.4 


18 ± 17 



CD 



0.1 



0.01 



0.001 
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11.9i?'? , for S > 5 mjy; 



• all z 
□ z < 1.7 
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0.01 0.1 
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Figure 7. Projected two-point correlation function for 100-fim- 
selected sources in the EGS field with Siooum ^ 5 mjy and z < 
1.7 (blue, open squares) compared to that produced by all sources 
in the EGS subject to the same flux cut (red, filled dots). The 
short-long dashed line shows the best power-law fit obtained for 
the more local sample by assuming 7 = 1.8 (see text and Table 2 
for details) 



n A (S > 5 mjy) ~ 0.91, n B (S > 5 mjy) ~ 0.09; n A {S > 
2 mJy) ~ 0.9, n B (S > 2 mjy) ~ 0.1. 

Once again assuming a power-law form with 7 = 1.8 for 
the spatial correlation function of low-redshift sources, the 
above values lead to clustering lengths respectively of: 
r (z < 1.7) = 4.lt?.o, for S ^ 8 mjy (only slightly changed 
from the previous analysis); 
r (z < 1.7) = 5.OI3 2 , for S ^ 5 mjy; 
r [z < 1.7) = 5.1±2:i. for 5 ^ 2 mjy; 
and minimum halo masses of: 
Logio[M min /M ]= ll.lig;*, for S > 8 mjy; 



-l.i) 
=12.0+ - 4 



Logi O [M min /M ] 

Logio[M min /M Q ]=12.0t^ 7 , for S > 2 mjy 

All the above values are summarized in Table 2. 

In all cases, although somewhat masked by the large 
uncertainties, the net result is a decrement of both the cor- 
relation length and of the halo mass of the host. Such a 
decrement is increasingly more pronounced as the fraction 
of high-z FIR galaxies increases. 

In order to cross-check the reliability of the above re- 
sults, in the EGS field we have also directly examined the 



clustering properties of the 592 Sioop 



5 mjy galax- 



ies endowed with either a spectroscopic or a photometric 
redshift determination, z < 1.7. This is presented in Fig- 
ure V7\ by the (blue) open squares and compared with the 
filled (red) points which reproduce the w(8) obtained for 
EGS galaxies at all redshifts. The amplitude decrement, al- 
ready appreciable by eye, implies in this case a correlation 
length ro = 5-llYg and a minimum mass for the host halo 
Logio[M m j n /M ] = 11.8i" l j , 'g, in excellent agreement with the 
previous determinations. 

A similar approach performed on z < 1.7 GOODS-S galax- 
ies with Sioo/jm ^ 2 mjy leads to ro = 4.3Jr 2 3 an d 
Mmin = II.5I3Q [M©]- Despite the uncertainties, the val- 
ues just found are once again in excellent agreement with 
those derived by making use of equations [5] and [6] 



5.2 Results and comparisons with other works 

The values obtained in Section 5.1 and Table 2 indicate that 
there is a substantial constancy in the clustering strength of 
relatively local, z < 1.7, sources as a function of their 100/im 
flux limit. 

We can actually extend the above analysis to create a 
homogeneous sample of galaxies selected at the same rest- 
frame wavelength in order to study the clustering evolution 
of star-forming galaxies free from the usual biases connected 
with the selection band. In fact, as it is possible to appre- 
ciate from the first five lines of Table 2, the lOO^im-selected 
samples at the various flux cuts and for 2 < 1.7 identify 
a redshift range (where the limits correspond to the rms 
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of the distributions presented in Figure Hb, which in turn 
corresponds to a range of rest-frame wavelengths all cen- 
tred at A ~ 60 fmx. It is then possible to compare such 
results with those obtained in the nearby universe for IRAS 
galaxies. This is done in Figure 151 which reports the red- 
shift evolution of the clustering strength of 60/xm-selected 
galaxies both by looking at their correlation length (left- 
hand panel) and minimum halo mass (right-hand panel). 
The filled points indicate the results of our work, while the 
open circle in the left-hand panel is the result for ro obtained 
by Saunders et al. (1992), once the slope 7 was changed from 
their measured value to the 7 = 1.8 used throughout this 
work. The corresponding point in the right-hand panel of 
Figure [8] for the minimum halo mass of IRAS galaxies was 
estimated with our code, to reproduce the correlation length 
reported by Saunders et al. (1992). The high, z ~ 2, redshift 
points, illustrated in Figure[8]by the open squares, represent 
the values for the correlation length and M m i n obtained by 
Magliocchetti et al. (2011) for a sample of Herschel-PEP 
galaxies, selected in the redshift range z — [1.7 — 2.6] at 
160/^m in the GOODS-S field. In fact, in the considered red- 
shift range, an observed A = 160^m corresponds to a rest- 
frame value of 52 ± 5/um (cfr Table 2), perfectly compatible 
with our chosen sample. 

What we have at hand is therefore a sample which in- 
cludes star-forming galaxies all selected at the same rest- 
frame frequency. Furthermore, most of the sources belong- 
ing to this sample have comparable bolometric luminosities 
and therefore star formation rates (SFR henceforth). In fact, 
we can use the redshift distributions from Section 3.2 to 
estimate the mean minimum bolometric luminosity which 
corresponds to the various flux cuts Sw&n via 



< L n 



with 



>= 



f*2Z L 'mU^)N{z) dz 

"7;°" n(z) dz ' 



L' min (z) = / / m m(A*,z) /(A) d\, 



(7) 



(8) 



where Z mm (A*,z) = ^TvS- m i- n d 2 L /K{\* ,z), is the minimum 
monochromatic luminosity at A* corresponding to Smin, d.L 
is the luminosity distance and the K-correction is expressed 
as K(\*,z) = (l+z)f(\*/(l+z))/f(\*). f is the normalized 
emission spectrum of the source. This was chosen in agree- 
ment with the results obtained by Gruppioni et al. (2010) 
for the population mixture of PEP-selected galaxies, to be 
of the kind generated by a moderate star-forming source 
(M82-like spectrum) at low-to-intermediate redshifts and as 
produced by an intense star-forming activity (Arp220-like 
spectrum) at z ~ 2. However, we stress that the following 
results do not greatly vary if we assume the same spectrum 
for all objects. 

The results of the above analysis are summarized in 
Table 2, which clearly shows that the values for the av- 
eraged minimum bolometric luminosities obtained at rela- 
tively low redshifts at the various flux cuts coincide within 
the errors, and, in addition, that the luminosity of 160/im- 
selected sources in GOODS-S at z ~ 2 is comparable 
with those derived for 100/im-selected, z < 1.7 sources in 
the COSMOS and EGS fields. By converting such lumi- 
nosities into star formation rates with the usual formula 



SFR[M Q /yr] = 1.8 x lO~ lo L bo i/L (Kennicutt et al. 1998), 
we then obtain that the high-z sample includes sources which 
are as active as low-redshift COSMOS and EGS galaxies, 
i.e. the sources included in our analysis are not only selected 
at ~ the same rest-frame wavelength, but also show compa- 
rable luminosities and star-formation rates. In other words, 
the sample of PEP-selected galaxies presented in this part of 
the work is reasonably free from both selection and luminos- 
ity biases, so that the evolution of their clustering properties 
can be considered as only driven by the intrinsic evolution 
of the galaxy population. / 

We note that the IRAS sample has been included in our 
analysis even though its sources are in general fainter than 
those observed by Herschel (even though with luminosities 
which are comparable to those of the z < 1.7 GOODS-S 
sample). This is because it has been shown (e.g. Hawkins 
et al. 2001) that their clustering properties do not depend 
on their luminosity (while they do depend on their colours), 
even in the most extreme cases. The Hawkins et al. (2001) 
conclusions are also confirmed by our analysis which shows, 
within the errors, a remarkable consistency of the clustering 
properties of 2 < 1 dusty star-forming galaxies with SFR. 

Figure 8 then shows the redshift evolution of the clus- 
tering properties of sources with comparable absolute lu- 
minosities and star formation rates and which are selected 
at roughly the same wavelength at the various cosmological 
distances. There is a clear trend for the clustering strength 
of these galaxies to increase at earlier epochs. This implies 
that also the masses of the hosts tend to increase as one 
moves to higher redshifts, showing a very clear evidence for 
downsizing. However, while the rise of both the correlation 
length and the halo mass is relatively limited up to z ~ 1, so 
much that it could be considered as almost constant within 
the uncertainties, there is a jump of about a factor 3 in ro 
and of a couple of orders of magnitude in the halo mass 
when moving from z ~ 1 to z ~ 2. The best linear fit to 
the data presented in the left-hand panel of Figure 8 is pro- 
vided by a functional form Logio[M m i n /Mo] = az + ft, with 
a - 1.1 ± 0.2 and (3 - 11.3 ± 0.2 and it is illustrated by the 
dashed line. 

Also, while the halo masses of FIR-selected galaxies in 
the relatively local universe all lie below MoPj (represented in 
the right-hand panel of Figure 8 by the dotted line), mean- 
ing they are relatively common objects, the halo mass found 
at z ~ 2 (M m i n as large as 10 4 Mo-) corresponds to extremely 
high density (and therefore rare) peaks of the underlying 
density field. Furthermore, Magliocchetti et al. (2011), by 
estimating the number density of dark matter halos via 
straightforward integration of the Sheth & Tormen (1999) 
mass function found that "every dark matter halo which 
reaches the mass to host a (dusty) star-forming galaxy at 
redshift z ~ 2 will be inhabited by (at least) one of such 
objects". 

This, and the results obtained so far in this work, im- 
ply that intense star-forming activity within galaxies takes 
place in very different environments at the different cos- 
mological epochs. At relatively low, z < 1.7 redshifts, the 
hosts of such systems are low-mass and isolated galaxies, 



1 M CT is the typical mass-scale at which a lcr fluctuation is col- 
lapsing. 
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Figure 8. Evolution of the clustering length ro (left-hand panel) and of the minimum halo mass M m i n (right-hand panel) for a sample of 
galaxies selected at about the same, A ~ 60/im, rest-frame wavelength. The filled points show the results obtained in this work, while the 
empty squares are taken from Magliocchetti et al. (2011) and are derived for a sample of z ~ 2 GOODS-S galaxies selected at 160/xm. 
The open dots are the Saunders et al. (1992) result for IRAS galaxies. In the right-hand panel, the dashed line represents the best fit 
to the data obtained for a functional form Log[M m i n /M0] = a ■ z + /3, with a = 1.1 and /3 = 11.3, while the dotted line represents the 
redshift evolution of Mo- (see text for details). 



and only a fraction (which decreases as one approaches the 
more local universe) of the virialized halos of comparable 
mass will be actively forming stars. On the other hand, the 
very same amount of star formation found at higher redshifts 
happens in extremely massive/proto-cluster like structures, 
and is presumably triggered by interactions of close pairs. 
Furthermore, this process is very common, as approximately 
all the z ~ 2 halos which reach the threshold mass to host a 
star-forming galaxy will be inhabited by at least one of such 
objects. 



6 CONCLUSIONS 

We have used data taken by the Herschel-PEP survey at 
100/mi in the COSMOS and Extended Groth Strip fields to 
investigate the possible dependence of the clustering prop- 
erties of FIR sources on their fluxes. To this aim, we have 
selected two samples of galaxies respectively brighter than 
8 mjy (COSMOS) and 5 mjy (ECS) and estimated their 
projected angular two-point correlation function w{0). 

By fitting the available observations to the functional 
form w[ff) = A6 X ~"' (with 7 fixed to the value 7 = 1.8) 
and by then de-projecting the best-fit values by means of 
the self-consistent redshift distributions in the two fields, we 
find for the correlation length: ro(Sioo^m ^ 8 mjy) =4.3JIq'7 
Mpc for COSMOS sources and ro(Siooi*m > 5 mjy) =5.8±2;o 
Mpc for ECS galaxies. These data show a mild trend for the 
clustering strength of FIR-selected sources to decrease with 
increasing fluxes. Such a trend is confirmed if we add to the 
above results those obtained by Magliocchetti et al. (2011) 
for PEP galaxies brighter than 2 mjy in the GOODS-South 



A similar (inverse) dependence on the flux limit is found 
for the minimum mass of the dark-matter halos hosting 
such FIR-selected sources. By applying the extended Press 
& Schechter theory (Sheth & Tormen 1999), we derive from 
our data: Logio[M min /M ](SiooMm ^ 8 mjy) 
Logio[M min /M ](S'iooMm > 5 mJy)=12.4 



1+0.4 

-0.6- 



11.6± ;2 and 



field MSioo/xm > 2 mjy) =6.3^3 Mpc). 



We also report a remarkable flattening of the observed 
w(9) for angular distances 0.004 < 6 < 0.01 degrees. Such 
a trend, which proceeds down to the smallest probed 
scales in the case of sources with fluxes Sioojun 1? 8 mjy, 
is due to the well known "avoidance behaviour" of local 
star-forming galaxies which prefer to stand isolated and 
not to have companions within the same dark matter halos 
(cfr Magliocchetti & Porciani (2003) or Zehavi et al. (2011)). 

But what is the physical explanation for the above findings? 
A careful investigation of the redshift distributions of PEP 
sources selected at lOOpim shows that the effect of lowering 
the flux limit is basically two-fold: on the one hand it moves 
the peak of the distribution to (slightly) higher redshifts, 
while more importantly it affects the appearance and rela- 
tive importance of a population of z ~ 2 sources. Such a pop- 
ulation makes about 10% of the total counts for Sioo/jm ^ 2 
mjy, while it already becomes negligible for sources brighter 
than 8 mjy (~ 4%). 

Magliocchetti et al. (2011) have shown this z ~ 2 
population to be very strongly clustered, with correlation 
lengths as high as ro ~ 19 Mpc and corresponding halo 
masses M m i n ~ 10 13 ' 9 Mq, so that it is plausible to assume 
these objects responsible for the enhancement of the 
clustering strength of faint, Sioo^m < 8 mjy FIR sources. 
We have then sub-divided our samples into intermediate 
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[z < 1.7) redshift galaxies and higher-z (z ~ [1.7 — 2.6]) 
ones and applied two different methods to evaluate the 
clustering properties of intermediate redshift FIR galaxies: 
one which relied on the subtraction of the clustering signal 
produced by the high redshift population from the total 
measured angular correlation function, and one which 
directly estimated w(&) for z < 1.7, 100/im -selected 
galaxies in the GOODS-S and EGS fields. Both methods 
converge, indicating an overall constancy of the clustering 
strength of relatively local sources, independent of flux cuts. 

With the above results in hand, we have then studied the in- 
trinsic redshift evolution of the clustering properties of dusty 
star-forming galaxies. This was done by extending our anal- 
ysis to comprise sources in the redshift range [0-2.5] which: 
a) were all selected at the same rest-frame wavelength to 
minimize selection biases; b) showed comparable luminosi- 
ties and therefore star-formation rates in order to overcome 
Malmquist-like bias effects. 

To achieve this, we added to our set of 2 < 1.7 100/im- 
selected sources both data from the clustering analysis of 
z ~ 0, 60/im-selected galaxies from the IRAS survey (Saun- 
ders et al. 1992) and also the results from Magliocchetti et al. 
(2011) for z ~ 2, 160/im selected galaxies in the GOODS-S. 
All these objects were observed at approximately the same 
rest-frame A = 60/im. Furthermore, the star formation rates 
estimated for z ~ 2 sources (SFR> 100 M0/yr) are com- 
parable with those derived for < z >~ 0.5 galaxies in the 
COSMOS and EGS fields. 

Our analysis for this homogeneously-selected sample 
shows a strong evolution of the clustering properties with 
look-back time. More specifically, both the clustering length 
and the minimum halo mass can be considered as approx- 
imately constant (or at most very weakly dependent on 
redshift), irrespective of the chosen flux cut, up to about 
z ~ 1. Above this redshift value, there is a huge jump of 
a factor ~ 3 in ro and of about two orders of magnitude 
in M m j n . In this latter case, the dependence of M m i n on z 
can be parametrized as Logio[M m i n /M0] = az + j3, with 



1.1 



+0.2 



and ft 



II 3 +0 ' 2 



This and the other results obtained in this work im- 
ply that intense star-formation within galaxies takes place 
in very different environments at the different epochs. In 
the relatively local universe, the hosts of such systems are 
low-mass and isolated galaxies, with only a fraction of them 
actively forming stars. On the other hand, the very same 
amount of star formation found at higher redshifts hap- 
pens in extremely massive/proto-cluster-like structures, and 
is presumably triggered by interactions of close pairs (cfr 
Kartaltepe et al. 2012). Furthermore, this process is very 
common, as all the z ~ 2 sources which reach threshold 
halo masses to host a star-forming galaxy will be inhabited 
by at least one of such objects. 

We note that such extreme findings for the clustering 
properties of star-forming galaxies at redshifts z ~ 2 — 3 
are not just the result of a particularly biased field or as 
a particular feature of the FIR realm. Indeed clustering 
length values which are comparable with those presented 
in this work and in Magliocchetti et al. (2011) are found 
in the same redshift range for high-SFR objects in various 
deep fields and at all probed wavelengths: from the UV (e.g. 
Magliocchetti & Maddox 1999), to the sub-mm regime (e.g. 



Hickox et al. 2012), passing through the Spitzer-2A[im. (e.g. 
Magliocchetti et al. 2008; Brodwin et al. 2008; Starikova 
et al. 2012) and the BzK selection (e.g. Lin et al. 2012). 
We will investigate the theoretical implications of such a 
behaviour in a forthcoming paper. 
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